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1 Introduction 

For weakly or moderately strongly correlated systems ab initio methods, such as the density 
functional formalism 01|2]| or the GW method JUS, are often quite successful. For strongly 
correlated systems, however, these methods are often not sufficient. It is then necessary to treat 
correlation effects in a more accurate way. Such systems are often quite complicated with large 
unit cells. It is then very hard to treat correlation effects within an ab initio approach, and one 
often turns to model Hamiltonians. The idea is then to focus on states and interactions believed 
to be particularly important for the physics of interest. This has the additional advantage that 
it may then be easier to understand the physics, since less important effects do not confuse 
the interpretation. On the other hand, there is a risk of oversimplifying the model and thereby 
missing the correct physics. The purpose of this paper is to discuss this approach. 
In principle it is straightforward to construct a model. We can produce a complete basis set and 
then calculate matrix elements of the real space Hamiltonian 



For atoms or small molecules, this Hamiltonian may then be solved using various many-body 
methods, e.g., configuration interaction (CI), where the many-body wave function is written as 
a linear combination of determinants. For strongly correlated solids, however, a Hamiltonian 
obtained in this way is often too complicated to allow reasonably accurate calculations. We are 
then forced to use substantially simpler models. This usually involves a drastic reduction of 
the basis set and the neglect of many interactions. Typical examples are the Anderson Q, the 
Hubbard [6] and the t - J models. 

This approach involves the neglect of interactions which are large. For instance, the Anderson 
impurity model is often used for a 3d impurity in a weakly correlated host. We define a direct 
Coulomb integral 



where ^j(r) is the wave function of a state i. Then the Coulomb integral F 3d 3d between 3d 
electrons is kept, while, for instance, the integral F^is between a 3d and a 4s electron is 
neglected. For a free Mn atom F 3 d,3d = 21 eV and F^is = 10 eV. Such an approximation 
is clearly highly questionable. An essential task is then to try to include explicitly neglected 
interactions or states implicitly as a renormalization of parameters in the model. As we show 
later, this leads to an effective Coulomb interaction between the 3c? electrons which is much 
smaller than the calculated value for a free atom. A basic assumption of such simple models is 
then that all the neglected interactions can with a reasonable accuracy be included implicitly as 
a renormalization of various model parameters. In this approach it is important to keep track of 
what effects are explicitly included in the model. These should not be included in the calculation 
of parameters, since this would involve double-counting. 

There are various ways of obtaining parameters. One approach has been indicated above. We 
use ab initio calculations to calculate parameters and then we try to estimate how these are 
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renormalized by neglected interactions. Another is to calculate certain properties of the model, 
compare with experiment and then adjust parameters until the experimental value is obtained. 
This approach then automatically gives renormalized parameters. It is important to try to obtain 
as much independent information as possible about the parameters, both from calculations and 
from different experiments, and to check if various pieces of information are consistent. 

The importance of obtaining theoretical information about parameters can be illustrated by the 
historical development of the theory of Ce compounds. Traditionally, Ce compounds were de- 
scribed in the so called promotional model [H. It was assumed that the Ce 4/ level was located 
very close to the Fermi level, E F , and that it had a very weak interaction with other states. A 
mean-field theory was then used to show that this leads to a very narrow resonance, as indicated 
in Fig. [Q The narrowness of the resonance could explain the large susceptibility and specific 
heat of Ce compounds, and the closeness of the 4/ level to Ep the change of apparent valence 
when the pressure or temperature were changed. Thermodynamic considerations, however, 
showed that the 4/ level ought to be about 2 eV below E F [9], in strong disagreement with 
the model. This result was later reconciled with experiment in a many -body approach [fTOlfTTIl . 
showing that even if the 4/ level is far below Ep it can form a Kondo-like many-body resonance 
at Ep leading to very large values of the susceptibility and the specific heat. This illustrates how 
an oversimplified (mean-field) method can nevertheless lead to reasonable results if it is com- 
bined with a bad choice of parameters. Correcting the parameters then forces us to use a better 
method and to find out more about the correct physics. 




Fig. 1: Schematic density of state for a Ce compound according to the promotional model. 
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2 Projecting out states 
2.1 One-particle Hamiltonian 

One approach to the construction of models is to project out states which are believed to to not 
be essential for the physics. We can illustrate this for a one-particle Hamiltonian 

# = E + E *0'^<Vj (3) 

i i^j 

We introduce a projection operator 

P = £>>H (4) 
where \v) are states we want to keep. We introduce the resolvent operator 

(* " H)- 1 = J2 W)("\(* - H)- 1 E MM = E I^TTF-H. ( 5 ) 

which has poles for z = ^ at the eigenvalues. Introducing the complement Q — (1 — P), we 
can write the Hamiltonian as H12U13I 



Hqp Hqq 



(6) 



where, e.g., H PP = PHP. Then we can derive the exact result 

P(z - H)' l P =[z- H PP - H PQ (z - Hqq^HqpY 1 . (7) 

The operator P(z — H)~ l P has the same poles as the original operator (z — H)^ 1 , if the 
corresponding eigenstates have weight inside the space P. The new operator has a smaller 
dimension, but because of the z dependence it is not simpler. To simplify the expression, we 
put z equal to an energy (e ) in the range of interest. The operator is then energy independent. 
As an additional simplification, we may assume that the off-diagonal elements of Hqq can be 
neglected. Then the matrix elements of the new operator become 

Uj -> Uj - E i lfl % ■ ^> 

This latter approximation is accurate if the states being projected out are much higher in en- 
ergy that the states of interest and if the off-diagonal elements are small compared with the 
energy difference e — E^. The assumption about Hqq being diagonal can also be relaxed. 
This approach reduces the size of the Hamiltonian matrix (reduces the number of states) at the 
cost of obtaining more long-range hopping. For a one-particle Hamiltonian, this approach is a 
controlled and systematic procedure for reducing the size of the Hamiltonian. 
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2.2 Many-body Hamiltonian 

We now consider a many-body Hamiltonian, with a two-body interaction in the form of a 
Coulomb interaction. We then define P as projecting out states that have no electron in cer- 
tain (high-lying) one-particle states and Q = 1 — P. We consider a Coulomb interaction 
with four (creation and annihilation) operators and project out a state with one electron in \/f). 
Then Plqp contains an operator c< an Hpq and operator c M(T . Even if we assume Hqq to be 
diagonal, we are left with an operator H P qHq P acting on a state without electrons in Then 
c Mo- c t CT = 1> an d two operators drop out. But we are still left with six other operators, which 
in the general case are all different. We have then generated a three-body operator. This is too 
complicated, and all such operators need to be neglected. Unless it can be shown that these 
terms are small, this means that there is not a controlled systematic procedure for reducing the 
number of states. We then have to rely on more intuitive approaches. 

As a simple example we consider a very simple model which is relevant for 3d impurities. The 
model is constructed so that an exact solution can be found. We want to illustrate how this 
model can be projected down to a simpler model with renormalized parameters. We introduce 
the Hamiltonian [14 J 



H = £iUia + (*^l<T^2ff + V ^L^4a + H- C -)] + U dd n 2 ^n 2i + U sd ^ U 2a n ia' ( 9 ) 

where level 2 corresponds to a 3d level and level 4 to a 4s level on a transition metal atom. 
Level 1 and 3 correspond to a ligand coupling to the 3d atom via the hopping integrals t and 
V. On the 3d atom there is a large Coulomb interaction Udd between electrons in the 3d level 
and a weaker U s d interaction between the 3d and 4s levels. We assume that orbital 2 is quite 
localized, so that t is small, but that levels 3 and 4 are delocalized, so that V is large. The level 
structure is shown schematically in Fig. [2j 



Fig. 2: Schematic picture of a very simple model of a transition metal compound, with a 3d 
atom (levels 2 and 4) coupling to a ligand (with levels 1 and 3). 

We first consider the spinless case, and put one electron in each of the spaces 1+2 and 3+4. 
We derive parameters in an intuitive approach, and then compare with a controlled projection 
approach, possible in this case. We introduce the eigenstates of the space 3+4 with the electron 
in space 1+2 on site 1 or 2. With the electron on level 1 the bonding and antibonding eigenstates 



i 



v 
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are 



Ipbl = ^3^3 + 

Ipal = ~ ^3^4, 

with the energies £ bl and e a \. With the electron in 1+2 on level 2 the states are 

^62 = (cos0)V>m + (sin0)^ a i 

^2 = (Sin0)-061 - (COS0)^ O 1, 



(10) 



(11) 
(12) 

(13) 



with the energies Eb2 and £ a 2- Here <p is of the order U s d/V which we is small in the limit 
we consider below. We assume that the electron in the space 3+4 can adjust completely to the 
movement of the electron in space 1+2 due to \V\ >> \t\. We then replace the four-level model 
in Eq. © by a two-level model with the effective level positions 



_eff 



eff 

fc 2 



^2 + £b2 



(14) 



To test this, we now solve the full model exactly. We introduce a complete basis set 

|i>=VWi|o> 



|2> 
13} 
14} 



v4vi|o> 



(15) 



where we have chosen the basis set so that only the first two states are relevant if the assumptions 
above are correct. We now calculate the resolvent operator 031 



[z-H) 



-1 



/ 



V 



- E\ — £&i 
—t COS0 



t sind> 



—t COS0 
Z — 62 — Sb2 

—t sin0 






—t sin0 

Z — E\ — E a \ 
—t COS(b 



t sin0 \ 


—t COS0 

- £2 — £a2 ) 



(16) 



We now focus on the upper left 2x2 corner and use Lowdin folding [12] to project out the two 
high-lying states. For instance, the 11 element takes the form 



H\\ — £1 + £&i 



£1 - £ a i)sm 



Ei - E al ){z - £ 2 - £ a 2) ~ t 2 COS 2 (, 



(17) 



For simplicity, we put £ 1 = e 2 and assume that the term t 2 cos 2 in the denominator can be 
neglected. Putting z ~ £\ + e&i, we then find that the correction term in Eq. (fT7l) is of the order 
t(t/V)(U sd /V) 2 . If |V| > |t| and \V\ > U sd , it is indeed justified to neglect the correction 
term. We then find that the level positions difference, sf 1 — £lf , have corrections to zeroth order 
in (l/V), due to £&i and £52- These corrections are included in our intuitive approach above. 
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Table 1: Ground-state energy E , occupancy of level 2, n 2 , and susceptibility \ of the spin- 
degenerate model ((9]). The parameters are s x = £ 2 = £3 = £4 = 0, t = 1, Um = 4 and 
U sd = 2. 



V 


c cff -eft' 
2 1 


U cti 


E Q + 2V 
Renorm. Exact 


n 2 

Renorm. 


Exact 


X 

Renorm. 


Exact 


1.0 


1.17 


3.18 


-1.05 


-0.95 


0.380 


0.364 


0.314 


0.312 


1.5 


1.39 


3.21 


-0.97 


-0.90 


0.339 


0.326 


0.266 


0.262 


2.0 


1.53 


3.29 


-0.92 


-0.88 


0.317 


0.307 


0.240 


0.237 


3.0 


1.68 


3.44 


-0.87 


-0.85 


0.292 


0.287 


0.214 


0.213 


4.0 


1.75 


3.55 


-0.85 


-0.84 


0.280 


0.277 


0.202 


0.201 


6.0 


1.83 


3.68 


-0.83 


-0.82 


0.268 


0.267 


0.190 


0.190 


10.0 


1.90 


3.80 


-0.81 


-0.81 


0.259 


0.258 


0.181 


0.181 


20.0 


1.95 


3.90 


-0.80 


-0.80 


0.252 


0.252 


0.174 


0.174 



Then there is a second order correction to the hopping integral due to cos</>. This correction is 
due to the fact that the electron in the space 3+4 cannot completely follow the electron in space 
1+2 in the optimum way. This correction is usually neglected. 

We now turn to the same model with spin degeneracy. The exact solution can then be obtained 
from a 16 x 16 matrix. In this case the analytical calculation is to complicated to illustrate what 
happens, and we focus on a numerical calculation. We first calculate the energy E(n 2 ) of the 
3+4 space as a function of the occupancy of level 2. We then obtain 

ef = £! + £7(0) 

ef = e 2 + E(l) (18) 
U eS = £7(2) + £7(0) - 2£7(1) 

in analogy with the spinless case. We then calculate the ground-state energy, £7 , the occupancy 
of level 2, n 2 and the spin susceptibility \ — —d 2 E (H)/dH 2 , where the model couples to an 
external magnetic field via the term —H(n 2 -\ — n 2 ^. The results are shown in Tabled! We have 
added a contribution 2V to the total energy £7o, since there would have been a trivial contribution 
—2V if there had been no interaction between spaces 1+2 and 3+4. As expected, the agreement 
between the approximate (Renorm.) and exact results improve as \V\ is increased. However, 
the agreement is surprisingly good even for V — t. 

3 Effective Coulomb interaction 

The essential point in the model in the previous section is that we can distinguish between two 
types of electrons, "slow" electrons (space 1+2) "fast" electrons (space 3+4), in the following 
referred to as "localized" and "delocalized". The idea is that the delocalized electrons are 
assumed to adjust in an optimum way to the movements of the localized electrons. We can 
then estimate effective parameters in a similar way as in the previous section. For each system 
we then have to decide which electrons are localized and included explicitly in the model and 
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Table 2: "Slow" ("localized") and "fast" (" delocalized" ) electrons for 3d and Af compounds 



System 


Localized 


Delocalized 


Af compounds 




5d 


3d compounds 


3d 


As, Ap 



which are delocalized and only included implicitly as a renormalization of the parameters. This 
is illustrated in Tabled For Af compounds the Af DFT band width is about 1/10 of the 5d 
band width, and we may reasonably talk about two types of electrons. For 3d compounds this 
distinction is much less clear cut. 



3.1 "Perfect screening" 

We now focus on the calculation of an effective Coulomb integral U eS , as an essential model 
parameter. We apply the approach in the previous section to real system. For that reason, we 
need to know how the energy of the system varies with the occupancy of, e.g., a 3d or Af level 
[Eq. (fT8l) l. Herring fTTSll estimated these energies using atom data, assuming that any change 
in the number of localized electrons on an atom is compensated by the opposite change in the 
number of delocalized electrons on the same atom. For a 3d metal this can be written as 

U = E{3d n+l As°) + E{3d n ~ 1 As 2 ) - 2E{3d n As 1 ), (19) 

where E{3d n As m ) is the energy of an atom (ion) with n 3d electrons and m 4s electrons. In this 
approach is is assumed that the variation in the number of 3d electrons is perfectly screened by 
a change in the number of 4s electrons. We refer to this as "perfect screening". 
A similar method was used by Cox et al. lfT6l who studied transition metals and Herbst et al. ifTTl 
who studied the rare earths. They performed Hartree-Fock calculations for renormalized atoms 
with Wigner-Seitz boundary conditions. 

3.2 Constrained density functional formalism 

Dederichs et al. ifTSl calculated U using a constrained density functional formalism. The func- 
tional for a 3d compound is written as 

E[nU = F[n] + J d\V ext {v)n{v) + p{j d 3 rn(r) - N} + /4{ J d?r^ d {r) - n^}. (20) 

Here F[n] describes the kinetic and potential energy of the system, V ext (r) is an external poten- 
tial, /i is the chemical potential and fi\ d is Lagrange parameter. n\ d is the number of localized 
electrons on site i, referred to as the central site in the following. A stationary point of the 
energy functional to density variations is searched 

dF ■ 
= — + V ext (r)+p + pl d P* d , (21) 

where P 3d is a projection operator acting on the localized electrons. From this a new Kohn- 
Sham equation can be derived, where ji\ d enters as an additional nonlocal potential acting only 
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on the localized electrons. fi l 3d is varied until the prescribed number of 3d electrons is obtained. 
This requires a definition of localized electrons. For instance, in methods based on an expansion 
in spherical waves in a region around each nucleus a natural definition can be introduced. This 
way an effective U cS is calculated, using a formula equivalent to Eq. (fT8l) . 
In the approach above, U cS contains a change of the kinetic energy of the electrons included 
explicitly in the model. This contribution needs to be subtracted. Hybertson et al. |fT9l did this 
by considering the model Hamiltonian in which U cS will be used, e.g., 

h = J2 %Vi< + E ueSn nna, (22) 

ija i 

where tij are hopping integrals. This model is then solved in a constrained mean-field the- 
ory, to simulate the constrained density functional approach. The energy as a function of the 
constrained occupancies is calculated, and U eS is varied until the constrained DFT result is re- 
produced. We refer to this as cLDA. Cococcioni and Giroconcoli [20] used a similar approach. 
An alternative approach was used by McMahan et al. fl2TTl and by Gunnarsson et al. [|22l . They 
performed a band structure calculation with a large unit cell ET1 - I231 . Then the hopping integrals 
from the orbital with localized electrons is cut off for the central atom in the unit cell. Then 
the occupation of the orbital can be trivially varied without a variation of the kinetic energy for 
hopping in and out of the orbital, since this energy is zero. Double-counting is also explicitly 
avoided, contrary to claims elsewhere ll24l . This method is referred to as "cut off' LDA. In a 
different method, the hopping between the localized orbitals and the delocalized orbitals was 
cut on all sites, not only on the central site ll2~8l . 



3.3 Constrained RPA 

A different approach was taken by Aryasetiawan et al. [|25l . They calculated the Coulomb 
interaction using a constrained random phase (RPA) screening. In RPA the polarizability is 
written as 

occ unocc 

P(r,r' :,) = ^^ ^(r)^(r')^(r)^(r') (23) 

* 3 

x( \ I ), 

u — Ej + Si + i0 + oj + Ej — Si — i0 + 

where ^(r) and £j are one-particle eigenfunctions and eigenvalues. Calculating the effective 
Coulomb interaction by using this screening would be incorrect, since it would involve double- 
counting. The Hubbard model explicitly allows localized electrons to screen the interaction be- 
tween localized electrons, and the use of Eq. (1231) would then lead to double counting. Aryaseti- 
awan et al. 11251 therefore excluded contributions to Eq. (|23T) where both i and j stand for Bloch 
states containing mainly localized states. For a transition metal compound they then excluded 
states within an energy window where the states are mainly of 3d character and for a rare earth 
compound a window where the states are mainly of 4/ character. The definition of the energy 
window involves uncertainties [|25l . This method is referred to as cRPA. 



Model parameters 



5.11 



Table 3: Contribution to U for a free Mn atom with the configuration 3d 5A As°' 6A Ap ' 70 . This 
corresponds to the configuration for Mn in CdTe. 



Unrenormalized (F°) 


21 A eV 


Relaxation of 3d orbital 


-5.2 eV 


Relaxation of 4s, Ap orbitals 


-2.2 eV 


Relaxation core, XC effects 


-1.2 eV 


Atomic U 


12.8 eV 



Table 4: Contribution to U for a Mn impurity in CdTe. 



On-site relaxation 


15.4 eV 


Charge transfer from Mn 


-7.6 eV 


Charge transfer to n.n. ligand 


-0.4 eV 


Solid state U 


7.4 eV 



3.4 Screening and breathing 

The definition of U can be approximately rewritten as 

U = E(n + 1) + E(n - 1) - 2E(n) « -i, (24) 

on 

where E{n) is the energy of the system with n localized electrons and e is the energy eigenvalue 
of the localized orbital and n is the occupancy. If the system were not allowed to relax, this 
would lead to U = F, where F [Eq. ©J is the direct Coulomb integral of the orbital. In reality, 
the charge density relaxes and the corresponding change in the electrostatic potential acts back 
on the orbital eigenvalue, reducing the shift as n is varied and leading to a renormalized U. 
We can illustrate this for the case of a Mn impurity in CdTe ll2~6ll . First a free Mn atom is 
studied [Table 0. The spherical part F.® d 3d of the direct Coulomb integral is large, 21 eV. The 
main renormalizing process is a breathing of the 3d orbital, where the orbital expands as the 3c? 
occupancy is increased H26U29H . This reduces U by about 5 eV. Breathing of the 4s, Ap and core 
orbitals contribute less. The net result is a reduction of U from about 21 eV to about 13 eV. In 
the solid there are similar breathing effects, reducing U to about 15 eV (see Table©. However, 
now there is additional charge transfer from the surrounding to the Mn atom, reducing the U 
by almost 8 eV. Charge transfer to near neighbors (n.n.) plays a smaller role. The result is 
reduction of U to about 7 eV according to this calculation. 

The breathing effect can also be understood from Slater's rules [|30l . According to these rules, 
the effective nuclear charge for a 3d orbital is Z* = Z — 18 — 0.35(n 3 d — 1), where Z is the true 
nuclear charge and n 3d is the 3d occupancy. This illustrates how the effective nuclear charge is 
reduced and the orbital expands as n 3d is increased. According to Slater's rules, the occupancy 
of the 4s and Ap orbitals do not influence Z* for the 3d orbital. This then suggest that the 
charge transfer in the solid to 4s and Ap should not influence breathing very much. This is also 
supported by a comparison of Tables [3] and HI 
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Table 5: Results for U for Fe as an example of a 3d metal and Ce illustrating a Af metal. 



System 


cLDA 


"cut-off" 


cRPA 


"perfect screening" 


Exp 


Fe 


2.2||20l 


6J231 


4(121 


2.7 una 


2[|3Tl[32l 


Ce 


4.5 [20] 


6J23 


3.2-3.3 [25] 


5E3 


5-7 [35] 



3.5 Results 

We now consider results obtained using the methods above for 3d and 4/ metals. Specifically, 
we consider Fe and Ce as examples of 3d and 4/ metals. The results are shown in Table [5] 
"Perfect screening" provides a rather good estimate for both Fe and Ce. The "cut-off" method 
gives a substantially too large U for Fe. It was found 1(231 that only about half the screening 
charge is on the Fe atom, as one would expect from the energetics of the screening process 
l(23l . It is then to be expected that U is substantially larger than the "perfect screening" result. 
cLDA gives a very good result compared with experiment, and actually somewhat smaller than 
"perfect screening". It is not clear why this result is so different from the "cut-off" method, 
and it would be interesting to study the screening in cLDA. The U in cRPA is a bit too large. 
Interestingly, the "cut-off" method gives a good estimate of U towards the end of the 3d series, 
e.g. for the cuprates 1(341 . 

For Ce "perfect screening" provides a fairly accurate estimate of U. The "cut off" method 
gives only a slightly larger U, in good agreement with experiment. In this case it is found 
that the screening charge on Ce is approximately unity (l23l . so it is not surprising that there is 
rather good agreement with "perfect screening". cLDA gives a U slightly smaller than "perfect 
screening" and U in cRPA is substantially smaller. It is not clear why cRPA implies such an 
effective screening and gives a U that is only roughly half the experimental estimate. 

4 Neglected renormalizations 

In this section we discuss two renormalizations of parameters, which are usually neglected. 
The purpose is not argue that these effects should be included. This could be done, but it would 
result in more parameters and the results would probably be less transparent. The purpose is 
rather to illustrate that the parameters of effective models contain complicated renormalizations, 
and that ab initio estimates of such parameters may neglect several such effects. The purpose 
is also to show that if we insist on rather simple model, which is advocated here, the effective 
parameter may actually be different for different properties. 

4.1 Configuration dependence of hopping integrals 

We already discussed in Sec. 13.41 that there is a substantial breathing of the localized orbital 
when the occupancy is changed. This changes the hopping integral into this orbital. In the 
LMTO method ((36ll . used here, the hopping integral V is related to a potential parameter A, 

V 2 ~A^ S -{MC,s)} 2 , (25) 
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Table 6: Potential parameter A for different configurations of Mn, Ce and U in non-spin- 
polarized calculations. The localized orbital is 3d (Mn), 4/ (Ce) and 5/ (U), and we consider 
a core hole in the Is (Mn), 3d (Ce) andAf (U) orbital. The occupancy of the localized and core 
orbital is rii and n c , respectively. We introduce n°, which is 5 (Mn), 1 (Ce) and 3 (U) and n° c 
which is 2 (mn), 10 (Ce) and 14 (U). All energies are in Ry. 



ni 


n c 


Mn 


Ce 


U 


71° - 1 


< 


0.0051 


0.0008 


0.0072 




< 


0.0085 


0.0019 


0.0091 


n° + 1 


< 


0.0129 


0.0038 


0.0112 


"? 


n° c -l 


0.0040 


0.0005 


0.0053 


n° + 1 


n-10 


0.0067 


0.0011 


0.0069 



where <f>i(C, s) is the value of the localized orbital at the Wigner-Seitz radius s. The localized 
orbital with the angular momentum / is solved for an energy C, which gives the logarithmic 
derivative — / — 1 at s. The value of A is shown in Table[6]for a few metals with and without a 
core hole [|27l . The table illustrates the strong dependence of the hopping on the configuration 
used to calculate A. For instance, if we want to describe how a host electrons hop into a Ce 
atom, should we then use the initial configuration or the final configuration to calculate A or 
should we use an average? Table [6] shows that the difference could even be as much as a factor 
of two. 

To address this issue, we temporary introduce an impurity model with two orbitals 11271 

# = &(r,n?) (26) 
d 

where A is chosen so that <p\ is normalized. By forming linear combinations of 0° and <j)\, we 
can obtain an appropriate orbital for different occupancies, i.e., describing breathing. In, for 
instance, an Anderson impurity model we then introduce a term leading to transitions between 
these two orbitals 

U ^(4moi>om« + H.c.)(no + n x - n°), (27) 

ma 

where n, = 2~2ma n ima- If the occupancy of the two levels adds up to n°, the orbital 0° is 
appropriate and there is no mixing of the orbital <\>\. For any other occupancy transitions to 4>] 
are induced and the system has the freedom to adjust to the occupancy. For Mn in CdTe we 
find that U = 0.16 Ry. The energies of the two orbitals are quite different, e = —0.45 Ry 
and E\ = 1.68 Ry. The model then tends to have two sets of states, one set at s$ and one set at 
e%. We can then project out all high-lying states, having a substantial weight in <fij. The result 
is then that we recover the normal Anderson impurity model, with just one localized orbital. 
But in this process the hopping matrix elements are modified. Since the mixing matrix element 
U/(e% — £o) = 0.08 1, this approach should be rather accurate. 

We can then answer the question of how to calculate these elements. Let us consider a host 
electron hopping into a configuration with n\ localized electrons, resulting in a configuration 
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with rii + 1 localized electrons. The projection procedure then shows that the matrix element 
should be calculated using n\ + 1 electrons, i.e., the end configuration [571 . For ni = this 
is easy to understand. In the initial state there is no localized electron and the extent of the 
localized wave function then plays no role. It is then natural that it is the wave function in the 
final configuration that matters. In a similar way it is the initial configuration that matters when 
an electron hops out of the localized orbital. 

We then should be using different hopping integrals for different experiment. For Ce com- 
pounds, for instance, f° — /^-hopping is particularly important for valence photoemission, and 
we would use nj = 1 for calculating these hopping matrix elements. For inverse photoemis- 
sion, we are often interested in the relative weights of the f 1 and f 2 peaks. We then need to 
distinguish between the calculation of the ground state and the calculation of the final states, 
resulting from the inverse photoemission process. In the ground-state calculation the important 
matrix elements would be calculated for ni — 1 and in the final state for ni = 2. For core level 
spectroscopies we should in addition include a core hole for the calculation of matrix elements 
for the final states but not for the ground-state. 

As argued above, this would lead to a complicated model. It seems questionable if the possible 
additional gain in physics would justify such a complicated model with additional parameters. 
However, the example illustrates one source of uncertainty in models with configuration inde- 
pendent hopping parameters. It also illustrates how parameters can be different for different 
experiments. 

4.2 Many-body renormalization of hopping integrals 

In Sees. 12.21 andl3l we discussed how the effective level energies and Coulomb integrals can be 
obtained by letting delocalized electrons adjust to the movements of localized electrons. This 
approach, however, raises questions about other many-body effects. One issue is the Anderson 
orthogonality catastrophe 11371 . Consider the case when delocalized electrons interact with a 
(truly) localized electron via a Coulomb interaction. Let us then change the occupancy of the 
localized level by one and let |0) and 1 1) be the lowest states of the delocalized electrons in the 
presence of or 1 localized electrons, respectively. The (0|1) = for an infinite system [|37l . 
One might then think that the hopping integrals should be reduced by such effects. When 
a delocalized electron hops into a localized level, all the other electrons would adjust their 
wave functions to the new potential. Then one might expect that the overlap (0|1) = enters 
the effective hopping integral. This is, however, not the appropriate comparison. Anderson's 
orthogonality catastrophe refers to the case when the localized electron is removed from the 
system. Here it hops into or out of delocalized states. The appropriate comparison is then X-ray 
absorption (XAS) or X-ray emission (XES). In addition to the Anderson effect there is then an 
exciton like effect, transferring spectral weight towards the Fermi energy. For instance, the XES 
spectrum looks like 

S(u;)~(^^)<*e{u-u ), (28) 

Ul — CiJq 
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Table 7: Energy lowering AE and occupancy of the d level n d in the exact calculation ("Ex.") 
compared with results of calculations for the model rf29l) with U sd = 0. The unrenormalized d 
level position was used for "Unre." and the calculated renormalized position for "Ren." and 
"XAS". For "XAS" the effective hopping integral was renormalized [Eq. 001)/ and for "Fit" 
both the level position and the hopping were adjusted to obtain the best agreement with the 
exact results. The parameters are t — 1, B = 5 N = 17 and N c \ = 9. 

— rV d 





u sd 


Ex. 


Ren. 


Unre. 


Fit 


XAS 


Ex. 


Ren. 


Unre. 


Fit 


XAS 


_calc 
E d 


E d 


+ fit 
''eff 


-1.5 


1 


1.33 


1.28 


1.66 


1.33 


1.31 


0.89 


0.91 


0.94 


0.89 


0.89 


-1.09 


-1.09 


1.12 


-1.5 


2 


1.12 


1.02 


1.66 


1.12 


1.08 


0.82 


0.87 


0.94 


0.83 


0.81 


-0.79 


-0.81 


1.18 


-1.5 


3 


0.98 


0.83 


1.66 


0.99 


0.94 


0.76 


0.81 


0.94 


0.78 


0.74 


-0.57 


-0.64 


1.21 


-1.5 


5 


0.83 


0.62 


1.66 


0.88 


0.78 


0.66 


0.70 


0.94 


0.69 


0.62 


-0.29 


-0.41 


1.30 


-1.0 


3 


0.64 


0.48 


1.20 


0.69 


0.62 


0.57 


0.55 


0.90 


0.55 


0.53 


-0.07 


-0.09 


1.31 


-0.5 


3 


0.42 


0.29 


0.78 


0.44 


0.41 


0.33 


0.24 


0.79 


0.31 


0.31 


0.43 


0.36 


1.22 


0.0 


3 


0.29 


0.21 


0.44 


0.30 


0.29 


0.18 


0.11 


0.50 


0.17 


0.17 


0.93 


0.76 


1.15 


10 


3 


.043 


.040 


.043 


.044 


.043 


.004 


.003 


.004 


.004 


.004 


10.9 


10.1 


1.00 



where a; is a typical energy and cu is the threshold energy. The exponent a is positive and 
determined by the phase shift due to the Coulomb interaction between localized and delocalized 
electrons. From Eq. (1281 we might then expect hopping integrals for states close to the Fermi 
energy to be enhanced. This would then in particular influence thermodynamic properties. 
To check these ideas we have considered the spinless model ll38l 

N N jj N N 

h = J2 e ^ + e ^ + 7^ + Hx -) + ir S E ^ n ^ (29) 

k=l * k=l k=l 1=1 

where we have introduced N delocalized states with the energies e k and a localized state with 
the energy e d . There is a hopping integral t connecting the localized and delocalized states. 
When the localized level is occupied the delocalized electrons feel a scattering potential U sd . 
The delocalized levels are equally spaced over an energy 2B. 

This model can be solved using exact diagonalization for finite N ll38l . We have used iV = 17 
and the number of electrons N e i = 9. Although this is far from an infinite system, Anderson's 
orthogonality catastrophe already has an effect. For B = 5, e d = —1.5 and U sd = 5, the overlap 
between the lowest states of delocalized electrons in the presence or absence of a localized 
electron is 0.85 < 1. We then calculate the energy lowering AE = E — e k , where E is 
the ground-state energy and the sum goes over the N d lowest states. We also calculate the 3c? 
occupancy, n d and the charge susceptibility Xc = —dn d /de d . The results are shown in Table |7] 
and 1 

We compare the exact results with several approximations [|38l . In all these calculations U sd 
was put to zero and its effects were approximately included via renormalized parameter. The 
column "Ren." shows results where e d was replaced by e^ alc = E (i) — E o (0). Here E (n d ) 
is the energy of the model as a function of the occupancy n d . In this calculation the hopping to 
the localized level was cut to avoid double-counting. The Table also shows results for unrenor- 
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Table 8: Same as for Table\7\but for calculating x- 

Xc 



e d 


u sd 


Ex. 


Ren. 


Unre. 


Fit 


XAS 


_calc 
E d 


E d 


j-fit 
r eff 


-1.5 


1 


0.12 


0.10 


0.05 


0.12 


0.13 


-1.09 


-1.09 


1.12 


-1.5 


2 


0.20 


0.19 


0.05 


0.20 


0.23 


-0.79 


-0.81 


1.18 


-1.5 


3 


0.27 


0.30 


0.05 


0.28 


0.32 


-0.57 


-0.64 


1.21 


-1.5 


5 


0.36 


0.55 


0.05 


0.38 


0.40 


-0.29 


-0.41 


1.30 


-1.0 


3 


0.47 


0.74 


0.12 


0.50 


0.47 


-0.07 


-0.09 


1.31 


-0.5 


3 


0.41 


0.41 


0.35 


0.43 


0.37 


0.43 


0.36 


1.22 


0.0 


3 


0.21 


0.14 


0.75 


0.22 


0.19 


0.93 


0.76 


1.15 


10 


3 


.0006 


.0005 


.0006 


.0006 


.0006 


10.9 


10.1 


1.00 



malized parameter ("Unre."). We then treated the and as fitting parameters, and adjusted 
these parameters to obtain the best possible agreement ("Fit") with the exact results. Finally we 
have performed calculations where the hopping matrix element to a level Ek 

[t eS (e k )} 2 = t 2 S(\e k -e F + u \), (30) 

was related to the X-ray absorption or emission spectra. [t c fj (£*:)] 2 summed over all states is 
unrenormalized, but the hopping parameters to states close to the Fermi energy, ep are enhanced 
at the cost of hopping to the band edges. In calculations with \t C R (sk)} we used the renormalized 
level position ef lc . 

We first compare the exact results with the unrenormalized and renormalized results. The renor- 
malization improves the agreement with the exact results substantially. For most parameter sets 
the agreement is relatively good. For U s d large and for \e^\ rather small, there are still substantial 
deviations. "XAS" shows the results when the hopping is renormalized as well, using Eq. ((30l) . 
There is then a substantial additional improvement, and the agreement is now generally a rather 
good. Finally, we have treated both the hopping and the level position as adjustable parameters. 
This gives only a marginal improvement and sometimes the results are even worse. This is 
remarkable, since the d-level position is now also a fit parameter and £^ alc is sometimes rather 
different from e^. On the other hand, hopping is energy-dependent, and "XAS" presumably 
describes this better than "Fit". This suggests that Eq. (l30l) gives a quite good renormalization. 

In the case of Ce the delocalized states are primarily of 5d character. According to the Friedel 
sum rule we can then estimate the phase shift as 5 ~ 7r/10. This then gives a singularity index 
of the order a ~ 0.1. For thermodynamic properties we may then expect an enhancement of the 
order of {Cj/T k ) ga , where T K is the Kondo temperature. For, e.g., CeCu 2 Si 2 T K = 0.001 eV 
and t 2 may then be enhanced by a factor of two, if we assume u ~ 1 eV. As discussed above, 
we do not advocate including these effects explicitly in a model. However, one should be aware 
that thermodynamic and spectroscopic properties may be renormalized differently. 
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5 Fullerenes 

In this section we discuss the parameters for a molecular solid. As an example we use fullerenes 
1591 . Similar work has been done for TTF-TCNQ flU. 



5.1 Hopping 
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Fig. 3: Levels of the Cqo molecule. The left-hand part shows the levels obtained by using a basis 
of one 2s and three 2p orbitals per carbon atom ( sp 3 ). The right-hand part shows the levels 
obtained by using just one radial 2p orbital per atom (2p r ). The numbers give the amount of 
radial 2p character (2p r ) in the full calculation (after Ref. 091/). 

The important levels in a Ceo molecule can be described in a tight-binding picture including a 2s 
and three 2p orbitals on each of the 60 C atoms. The corresponding molecular levels are shown 
in Fig. [3] The molecule forms approximate sp 2 hybrids on each C atom which point towards 
the neighboring C atoms and radial orbitals p r pointing out of the molecule. The former orbitals 
interact strongly and form bonding and antibonding molecular orbitals at the lower and upper 
end of the spectrum, respectively. The p r orbitals interact much less and form molecular orbitals 
in the middle of the spectrum. The figure illustrates that these orbitals can be described rather 
well using only the p r orbitals. In the neutral molecule all orbitals up to and including the h u 
orbital are filled. 

Ceo molecules condense to a solid of rather weakly bound molecules. Thus the distance (~ 3 
A) between the closest C atoms on two neighboring molecules is much larger than the distance 
(~ 1.4 A) between two C atoms on the same C 60 molecule. The molecular levels then essentially 
preserve their identity in the solid, but the discrete molecular levels are broadened to narrow 
essentially nonoverlapping band. The alkali-doped fullerenes are of particular interest. In these 
systems the t\ u band is partly filled. Therefore the three-fold degenerate ti u molecular level is 
particularly interesting. 
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The band structure can be described in a tight-binding (TB) scheme. We first form a molecu- 
lar orbital corresponding to the t\ u level. The hopping between the molecules is described by 
hopping integrals V ppa and V ppn corresponding to hopping between orbitals pointing directly 
towards each other or orbitals pointing perpendicular to the connecting line of the centers. Fol- 
lowing Harrison, PTI we assume that the ratio of the n— and a-integrals is -1/4. Then 

*W = ^#e-^-*°); ^ = -~ ^o = 3.1 A, (31) 

where R is the separation of the carbon atoms. The prefactor R has been included to simulate the 
r-dependence of a 2p orbital as described by Slater's rules [|30l . The overall hopping strength, 
determined by v a , is adjusted to the band width in a band structure calculation, and the decay 
length A is determined from the dependence of the band width on the lattice parameter. Here 
we use the parameters p2ll43l 

A = 1.98 A _1 and v a = 0.917 eV. (32) 

The resulting TB band structure is compared with an ab initio band structure calculation in 
Fig. SI The agreement is quite good. The resulting band structure has a simple parameter- 
ization [|42ll44l . The dominating hopping between two molecules in this structure is given by 
two equivalent hopping integrals, with all other hopping integrals being substantially smaller. 
Effectively, we have therefore adjusted this parameter requiring that the TB band width should 
agree with the LDA band width. The shape of the band structure in Fig. |4]is therefore primarily 
determined by the geometry of the C 60 molecule and by the relative positions and orientations 
of the C 60 molecules in the Fm3 symmetry. 



Fig. 4: Band structure for a Cqq solid in the Fm3 structure (a) according to an ab initio LDA 
and (b) according to a TB calculation (after Gunnarsson et al. I\43\l ). 

5.2 Coulomb interaction 

We first consider the Coulomb integral Uq between two ti u electrons for a free Ceo molecule. A 
very simple estimate is obtained by assuming that charge density of the t\ u orbital forms a thin 
shell of charge on a sphere with the radius R ~ 3.5 A. Then Uq = e 2 /R = 4 eV. This neglects 
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that the orbitals breath when the occupancy is changed. To obtain a better estimate one can 
calculate how the £ lu eigenvalue changes with occupancy, using Eq. (l24l) and LDA. This leads 
to values of the order 2.7-3.0 eV Il4"5l447ll . Uq for a free molecule can also be estimated from 
experimental results using 

Uq = I p {C m ) - A{C m ) = E (2) + E (0) - 2E (1), (33) 

where E (n) is the energy of a C 60 molecule with n ti u electrons. This leads to Uq ~ 2.7 
eV [|39ll48l . in fairly good agreement with theory. 

We next consider U for a C 60 solid, following Antropov et al. H71 . U is strongly screened by 
the polarization of the surrounding molecules. To describe this we put the C 60 molecules on an 
fee lattice and assign a polarizability a to each molecule. An electron is added to the central 
molecule, and the surrounding molecules are allowed polarize in a self-consistent way. This 
polarization acts back on the electron and reduces the energy increase of the t\ u level by an 
amount SU. The summation over neighboring molecules is extended until it is converged. The 
U for the solid is then 

U = U - SU (34) 

The value of a can be determined from the experimental value of the dielectric function (4.4). 
llBTI Using the Clausius-Mossotti relation and the lattice parameter a = 14.04 A, this leads to 
a = 90 A 3 . Ab initio calculations using the density functional formalism gave a = 83 A 3 ll46ll . 
Using a = 90 A 3 , Antropov et al. 03 found SU = 1.7 eV. Together with U = 2.7, this gives 
U =1.0 eV. These values of U do not include the metallic screening from the t% u electrons in 
A 3 C 6 o compounds, and they are appropriate for models where the metallic screening is treated 
explicitly when solving the corresponding model. 

We next consider the nearest neighbor interaction V, which is obtained by calculating the in- 
crease of the energy of a ti u orbital on a molecule 1 when an electron is added to a neighboring 
molecule 2. This leads to the result 

V = e 2 /R - 5V, (35) 

where R is the nearest neighbor separation and —5V is the lowering of the t\ u orbital on 
molecule 1 due to the polarization of the surrounding molecules when an electron is added 
to molecule 2. For a = 14.04 A, Antropov et al. 11471 estimated that <5V=1.12 eV, resulting in 
V = 0.3 for the polarizability a=90 A 3 . The same value V = 0.3 eV was also obtained by 
Pederson and Quong ll46l . We can see that U is indeed substantially larger than V, and that it is 
justified to focus on the effects of U at first. 

U can be estimated experimentally from Auger spectroscopy [|49l[50l . A carbon Is electron 
is emitted in a photoemission process. This is followed by an Auger process, where a carbon 
2p electron falls down into the Is hole and another 2p electron is emitted. For noninteracting 
electrons, the Auger spectrum is just the self-convolution of the photoemission spectrum. For 
the interacting system, the Auger spectrum is expected to be shifted due to the interaction of 
the two holes in the final state. Indeed, Lof et al. [|49l found good agreement with the self- 
convoluted curve when this was shifted by 1 .6 eV. The experimental estimate of the Coulomb 
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interaction is then U = 1.6 ± 0.2 eV [49 1 as an average over all orbitals and about 1.4 eV for 
the highest occupied orbital. Since Auger is rather surface sensitive, this number may be more 
representative for U at the surface. One can estimate that U at the surface is about 0.3 eV larger 
than in the bulk, due to fewer neighbors and less efficient screening [|47l . This suggests that the 
bulk value of U for the t\ u and h u orbitals may be on the order U — 1.1 eV, which close to the 
theoretical estimate. U has also been estimated for K 6 C 60 in a similar way [50J, and giving a 
similar value U — 1.5 eV. 

5.3 Electron-phonon interaction 




5 20 50 100 200 500 1000 2000 

FREQUENCY (cm ') 

Fig. 5: Schematic representation of various phonons in A3C60 compounds. The figure shows, 
from left to right, (a) librations, (b) intermolecular Cqq-Cqo phonons, (c) A-Cqq phonons and 
(d)-(e) intramolecular H g modes. The figure indicates the radial and tangential character of the 
low-lying and and high-lying H g modes, respectively (After Hebard /l52l/). 

The electron-phonon interaction plays an important role for many properties of alkali-doped ful- 
lerides. For instance, superconductivity is believed to be due the electron-phonon interaction. 
Fig. [5] indicates the different types of phonons in alkali-doped C 60 compounds. The low-lying 
modes are librations (4-5 meV) and intermolecular modes (energies up to about 17 meV) in- 
volving alkali-Ceo and C60-C60 modes. The high-lying modes (34-195 meV) are intramolecular 
modes, where the molecules are deformed. All the low-lying modes have a rather weak cou- 
pling to the ti u electrons, and the main coupling is to the intramolecular phonons. Here, we 
therefore focus on the the coupling to these phonons. These phonons couple primarily to the 
level energies in contrast to the intermolecular phonons which couple to the hopping integrals. 
The C 60 molecule has 60 x 3 — 6 = 174 intramolecular modes. For symmetry reasons, however, 
the t\ u electrons only couple to modes with A g or H 3 symmetry. There are eight five-fold 
degenerate H g modes and two nondegenerate A g modes. The coupling to the t\ u level takes the 
form (51 

8 5 3 3 

H d - ph = J2 J2 Yl m'V'LV'm'AM + & Lv/) 
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(36) 



u=9 



-I m =1 



where creates a t\ u electron with quantum number m and b\ M creates a phonon in mode v 
with quantum number M. The first eight modes are H g Jahn-Teller phonons and the next two 
A g phonons. The coupling constants are g u and the coupling to the H g phonons is given by the 
matrices 
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and the coupling to the A g phonons by 



V 



(i) 



/ 1 o \ 
1 
1 



(38) 



The corresponding dimensionless electron-phonon coupling constant is [53 



10 



a = -n(o) Yi^ + ~mo) y 

3 v ; ^ fiw„ 3 v ' ^ hu} tJ 



(39) 



u=9 



where N(0) is the density of states per spin and molecule and co u is the frequency of the mode 
v. 

The theoretical calculation of the electron-phonon coupling for a solid is very complicated. 
Lannoo et al. Il54l showed that for intramolecular modes in fullerides, important simplifications 
follow from the large difference between the intramolecular (Ej) and intermolecular (W) energy 
scales. The coupling for a solid can then be obtained approximately from a calculation for a 
free molecule and the density of states N(0) of the solid. Thus, it is sufficient to calculate the 
shift Ae va of the t\ u levels a for a free C 60 molecule per unit displacement of the z/th phonon 
coordinate. One then finds that 



A 



r 2 
' va 



(40) 



This gives a molecular specific quantity which is multiplied by iV(0). Table |9] shows results for 
the electron-phonon coupling. The theoretical calculations by Antropov et al [55 J, Faulhaber et 
al. 11561 and Manini et al. 11571 are based on ab initio LDA calculations. The work of Iwahara 
et al. is based on the B3LYP functional with some Hartree-Fock exchange mixed in. There are 
substantial deviations between the distribution of coupling strength to the different modes in 
the different calculations. This distribution is very sensitive to the precise form of the phonon 
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Table 9: Partial electron-phonon coupling constants X u /N(0) (in eV) according to different 
theoretical calculations and derived from photoemission and Raman scattering. The energies 
uo v (in cmT 1 ) of the modes for the undoped system are shown. 
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eigenvectors. The deviations between the total coupling strengths are smaller. The work of Iwa- 
hara et al. gives a stronger coupling than the other three calculations. This is not so surprising, 
since this work is based on a rather different functional. 

An experimental method for estimating the electron-phonon coupling is the use of photoemis- 
sion data. Because of the relatively strong electron-phonon coupling, we expect to see satel- 
lites due to the excitation of phonons. The weights of the satellites give information about the 
strength of the coupling. This is essentially the Franck-Condon effect, but because of the Jahn- 
Teller effect the calculation of the satellite structure is rather complicated. The photoemission 
spectra of K 3 C 60 and Rb 3 C 60 have been analyzed along these lines (621 . Due to the broadening 
effects in a solid and due to the complications in the theoretical treatment of bands with disper- 
sion, however, it was not possible to derive reliable, quantitative values for the electron-phonon 
coupling. 

Photoemission spectra have also been measured for free C^ molecules. In this case the the- 
oretical treatment is substantially simpler (581 . In these experiments, a beam of Cg ions was 
created and a photoemission experiment was performed using a laser light source (htu =4.025 
eV) and a time of flight spectrometer. The spectrum resulting from emission from the t\ u level 
was measured. To analyze the results, we use the coupling in Eq. ((361) of the t lu level to the 
two A g and the eight five-fold degenerate H g modes. For this model the ground-state can be 
calculated by numerical diagonalization to any desired accuracy [ 58l . Furthermore, within the 
sudden approximation (63l . the photoemission spectrum can easily be calculated. A set of cou- 
pling constants are then assumed and the resulting spectrum is compared with experiment. The 
coupling parameters are varied until good agreement with experiment is obtained, thereby pro- 
viding an estimate of the couplings. The resulting spectrum is compared with experiment in 
Fig. [6]and the corresponding parameters are shown in Table[9] An uncertainty in this approach 
is that with the available resolution, it is not possible to distinguish between the coupling to 
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A g modes and H g modes with similar energies. The couplings to the A g modes were therefore 
taken from a calculation ll55ll . With this assumption, the couplings to the H g modes can then 
be determined. An equally good fit can, however, be obtained using other couplings to the A g 
modes if the couplings to the H g modes are changed correspondingly. 



Fig. 6: Experimental (dots) and theoretical (full line) photoemission spectrum of Cg . The 
theoretical no loss (dashed), single loss (dotted) and double loss (dashed-dotted) curves are 
also shown. The contributions of the different modes to the single loss curve are given by bars 
(H g : open, A g : solid). The inset shows the experimental spectrum over a larger energy range 
(after Gunnarsson et al. /ESl/j- 

Substantially later the experiments in Ref. Il58l were repeated by Wang et al. [|59l . It was now 
possible to obtain a better resolution. These data have been analyzed in a similar way as in 
Ref. [58J by Iwahara et al. [60J. There results are shown in Table [9] The total coupling is 
weaker than in Ref. [58J, but still substantially larger than in the ab initio LDA calculations. 
The agreement with the calculation using the B3LYP is better. 

Raman scattering provides a different method of estimating the coupling strength. The electron- 
phonon coupling allows phonons to decay into an electron-hole pair in the metallic fullerides. 
This decay contributes to the width of the phonon and can be measured in Raman scattering. 
Other factors may also contribute to the width, but one can try to eliminate these by subtracting 
the width of the phonons for a nonmetallic system, where a decay in electron-hole pairs is not 
possible. This was done by Winter and Kuzmany lloTi . and Table |9] shows results adapted ll39l 
from the experiments ED- The total weight does not differ much from what was obtained 
from photoemission [60J , but the distribution of weight between the different modes differs 
dramatically. Theoretically, it is found that in the solid there is a transfer of weight to lower 
modes, due to the coupling to electron-hole pairs ll64ll . This mechanism is operative for the 
Raman data but not for the photoemission data (taken for a free molecule). This may explain 
some of the discrepancy between the PES and Raman data. 
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6 Conclusions 

For complicated systems with strong correlation effects it is often not possible to obtain accurate 
ab initio solutions, but it is instead useful to turn to models. An important issue is then how to 
obtain parameters and how to renormalize parameters to include as much physics as possible. 
We have discussed how the basic principle is to try to include implicitly as a renormalization 
of parameters all effects not explicitly included in the model. On the other hand, we should 
not allow effects included explicitly in the model to renormalize parameters. For many -body 
systems there is no general systematic and controlled way of doing this. The basic assumption is 
often that the electrons can be put into two groups of "fast" (delocalized) and "slow" (localized) 
electrons, where the "fast" electrons are assumed to adjust to the "slow" electrons, and therefore 
can projected out. Such a division is, however, often not very clear cut. Nevertheless some 
methods have been relatively successful in obtaining parameters for certain classes of systems. 
We have, however, shown simple examples of many-body effects that are usually not included, 
but can have an appreciable effect on the parameters. In particular, renormalization effects 
may work differently for different experiment. We have also argued that it is important to try 
to extract parameters from different sources, both theory and experiment, to obtain a better 
understanding of the accuracy of the parameters. 
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